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ABSTRACT 


A  Newton  method  using  interval  arithmetic  for  an 
n-dimensional  vector  equation  f(x)  =  0  is  treated.  Some  of 
the  properties  of  the  Newton  operator  are  investigated.  This 
leads  to  two  criteria  for  the  non-existence  and  for  the 
existence  of  a  zero  of  f (x)  in  a  given  vector  interval. 
Furthermore  a  posteriori  error  bounds  for  a  given  real  vector 
sequence  converging  to  a  zero  of  f(x)  can  be  established.  An 
example  deals  with  the  application  to  convex  programming. 


ON  THE  NEWTON  METHOD  IN  INTERVAL  ANALYSIS 
Karl  Nickel 


1.  Introduction. 

Let  f (x)  be  a  real  valued  n-dimensional  vector  function  of 
the  n-dimensional  real  vector  x.  The  solutions  £  of  the  problem 

f(x)  =  0  (1) 


are  treated.  The  Newton  method  for  (1)  defines  a  sequence 
{x*v^}  by  setting 


x(0>  "appropriate", 
x(v+1)  :  =  Nx(v)  for  v 


(2) 


where  the  Newton  operator  N  is  defined  by 

Nx  :=  x  -  f’"1(x)  f  (x)  (3) 

(f1  is  the  Frechet  derivative  (Jacobian  matrix)).  Under  certain 
assumptions  (St  exists,  f'”^  exists  in  a  neighbourhood  of  St,  x ^ 
is  "close"  to  x,  ...)  the  sequence  {x^}  converges  super- 
linear  ly  or  even  quadratically  to  x.  Unfortunately  this  is  not 
true  for  general  starting  vectors  x^,  as  the  following  one 
dimensional  example  shows : 

Example:  Let  0  <  a  <  b  and  define 
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r  3x_l  x2 

2  Si  "  2  ^i3 


for  0  £  x  <,  a  ; 


f  (x)  :=  / 


for  a  <  x  <,  b  ; 
3/T7S  b1/6  (x-b/3)ly/3  for  b  <  x 


-f  (-x) 


for 


x  <  0  . 


V. 


There  exists  exactly  one  solution  £  =  0  to  (1)  in  (-•*>,  »)  . 

The  function  f  e  C^-00,  <=°)  with  f*  >  0  and  f'-3"  >  0,  hence  the 
operator  N  from  (3)  can  always  be  applied  with  the  result 

0  <  |x(0)|  <  a:  |x(v+1)/x(v) |  <  1-3 (1- |x(0) |/a)/ (3-2 |x{0) |/a) 

<  1,  hence  {x^v3}  is  always 

(v) 


a  <  |x(0) |  <  b: 


convergent  with  x  ■+•  St  »  0  ; 
x<2v),x(0)  j[<2v+1)  .  .x(0),  fx<v)j  is  an 

alternating  sequence  with  |x^v3|  =  |x^| 
and  is  not  convergent ; 

_ _  |x(v+1,/x(v) |  >  l+(l-b/|x(0) |)  >  1  : 

{x^v3}  is  always  divergent  with 
[  x  ^ v  ^  |  -*■+“>  for  v  ■+  “  . 

R.  E.  Moore  [4]  first  suggested  replacing  the  real  Newton 
method  (2) ,  (3)  by  the  following  interval  Newton  methods 


b  <  x 


(0) 
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appropriate" 


x 


(0) 


k(0>  - 


*\ 


x(v+1)  :=  x(v)  n  Nf  (x(v),  x(v))  ,1  (4) 

x(v+D  *  x(v+1)  arbitrarily  , 

> 

where  the  interval  Newton  transformation  operator  is  defined  by 

NT(x,  x)  :=  x  -  f'“1(x)  f (x)  .  (5) 

Here  {x^}  denotes  a  n-dimensional  interval  vector  sequence, 

F'  (x)  is  an  interval  extension  of  the  Jacobian  matrix. 


The  following  surprising  fact  seems  not  to  be  widely  known: 
In  the  one-dimensional  case  n  =  1  the  interval  Newton  method 
defined  by  (4)  and  (5)  never  fails  1  It  delivers  always  two 
convergent  sequences:  x^  -*•  St;  x^  x  with  x,  x^  c  x^  for 
v  ■  0,  1,  ....  if  only  f  t  C^,  if  F'  1  exists  and  if  x  t  x^ 

(see  K.  Nickel  [5]  and  R.  Krawczyk  [3]).  Hence: 

a)  While  the  real  Newton  method  (2) ,  (3)  gives  convergence 
only  in  the  neighbourhood  of  x,  the  interval  Newton 
method  (4),  (5)  provides  convergence  for  all  initial 
values  x^,  if  only  x ^  is  chosen  such  that  x  €  x^  . 

b)  The  interval  Newton  method  (4) ,  (5)  provides  not  only 

an  approximating  real  sequence  {x^}  with  x^  -*•  x,  but 
computes  at  the  same  time  a  sequence  of  error  bounds 
{x(v)}  with  x(v)  -*■  x.  In  the  case  of  the  real  Newton 
method  (2),  (3)  no  such  error  bounds  are  generally 
available. 


*1136 


-3- 


Applied  to  the  above  example  this  shows  that  while  the 
real  Newton  method  (2),  (3)  failed  if  |x^|  >  a,  the  interval 
Newton  method  (4) ,  (5)  gives  a  convergent  approximating  sequence 
{x(v)  }  and  convergent  error  bounds  {x^},  if  only  0  =  x  e  x^  , 
regardless  of  the  initial  value  x^. 

Unfortunately  a  similar  statement  is  not  true  in  the  case 
of  more  than  one  variable:  if  n  >  1  then  in  general  (4),  (5) 
does  not  give  two  converging  sequences  {x^}  and  {x^},  as 
counterexamples  show. 

The  purpose  of  the  following  paper  is  twofold:  to 
investigate  some  of  the  properties  of  the  operator  NT  and  to  find 
a  way  to  guarantee  convergence  of  the  sequence  x^  x.  It  will 
be  shown  that  two  very  simple  criteriae  can  be  derived  from  the 
Newton  operator  which  -  when  fulfilled  -  guarantee  either  the 
existence  and  uniqueness  of  a  solution  x  of  (1)  in  a  given 
interval  vector  x  or  guarantee  the  nonexistence  of  such  a 
solution  x  in  x.  Both  criteria  can  easily  be  programmed  for  a 
computer.  In  order  to  get  unrestricted  convergence  of  {x^v  },  a 
new  approach  is  recommended:  Since  the  definition  of  the  real 
sequence  {x^}  by  (4)  does  not  always  lead  to  a  convergent 
series,  it  is  suggested  to  compute  this  real  sequence  {x^}  by 
any  method  providing  convergence  x^  x.  Then  the  sequence 
{x(v)}  constructed  a  posteriori  by  (4)  is  a  sequence  of  error 
bounds,  i.e.  x  «  x^  and  it  follows  easily  that  it  always 
converges,  i.e.  Him  x^  =  x.  This  leads  to  very  satisfying  a 
posteriori  error  bounds  in  linear  lgebra  problems,  for  eigen- 
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value  problems,  for  the  finding  of  roots  of  polynomials  etc.  As 
a  special  example  it  will  be  shown,  that  this  approach  can 
successfully  be  used  for  the  purpose  of  convex  programming.  A 
simple  example  will  show  how  it  works.  Since  the  definition  and 
use  of  the  interval  arithmetic  and  interval  analysis  is  not  yet 
known  widely  enough,  special  emphasis  has  been  given  to  listing 
all  assumptions  and  definitions. 

In  what  follows  a  slightly  different  interval  Newton 
operator  N(x,  x)  will  be  used,  which  is  a  generalization  of 
NT (x,  x)  from  (5).  This  goes  back  to  an  idea  of  E.  R.  Hansen  [2]. 

2.  Notations. 

Let  R  denote  the  field  of  real  numbers.  Let  a,  3  t  R 
with  a  <  6  and  denote  by  y  :=  [a,  8]  the  closed  interval  with 

endpoints  a  and  8.  Let  E  (y)  be  the  set  of  all  closed  intervals 
in  y  with  endpoints  in  R;  denote  analoguously  by  E(R)  the  set 
of  all  intervals  in  R  . 

Let  4,  n,  C  «  R  with  n  <  C  and  let  X  :■  tn,  U  c  I  (R). 

The  special  function  intval  (£,  i)  with  Interval : 

R  *  E  (R  )  -*•  I  (R  )  is  defined  by 

X  if  4  «  X  ? 

intval  (4,  T)  :■  <  In,  £]  if  4  >  C  ; 

[4,  Cl  if  4  <  n  , 

L 

i.e.  intval  (4,  X)  is  the  smallest  interval  in  E(R)  containing 


#1136 


-5- 


both  £  and  £. 


Let  R  n  be  the  n-dimensional  Euclidean  vector  space  with 
elements  a,  b,  . ..,  z  and  the  corresponding  components  ay,  b^,  . .., 
zv  e  R  for  v  =  l{l)n.  Interval  vectors  and  interval  matrices 
are  vectors  and  matrices  with  intervals  from  I  ( 1R)  as  components. 

In  what  follows  they  will  be  distinguished  from  real  vectors  and 
real  matrices  by  a  bar;  hence  x  denotes  a  real  vector,  x  an 
interval  vector;  X  resp.  X  is  a  real  resp.  interval  matrix.  An 
n  x  n  interval  matrix  A  is  said  to  be  invertible  and  to  have  an 
inverse  B,  iff  any  real  n  *  n  matrix  A  «  A  has  an  inverse  A  ^ 
and  if  A  t  B.  Kindly  note  that  an  interval  inverse  in  general 
is  not  unique.  Let  En(Rn)  be  the  set  of  all  interval  vectors 
a,  b,  ...,  2.  with  the  corresponding  components  ay,  b^,  ..., 
zv  «  E(R)  for  v  =  l(l)n.  Since  real  vectors  can  be  treated  as 
degenerate  interval  vectors  it  follows  that  Rn  C  In(Rn).  Let 
x  t  Rn,  x  (In(Rn)  with  components  xy  and  xy  for  v  =  l(l)n. 

The  special  function  Intval  (x,  x)  with  Intval  :  Rn  x  En(Rn) 

-*■  In(]Rn)is  the  n-dimensional  analogue  of  intval  (g,  f)  and  is 
defined  by  y  :=  Intval  (x,  x)  «  En(Rn)has  the  components 
yv  :=  Intval  (x^,  xy)  for  v  =  l(l)n.  The  interval  arithmetic 
operations  in  ln(Rn)  are  defined  as  usual  (see  R.  E.  Moore  [4]); 
they  will  be  not  distinguished  from  the  real  arithmetic 
operations  since  no  confusion  should  occur. 

Let  u  «  Sn(Rn)be  the  special  interval  vector  with  components 
uv  :=  y  :=  [a,  83  for  v  =  l(l)n.  Hn(u)  is  the  set  of  all 
interval  vectors  in  u  (i.e.  whose  components  are  contained  in  y) . 
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Let  f(x)  =  f(x^,  x2/  •••/  xn)  with  f:  u  -*■  R n  and  F(x)  = 

7(xv  3«2»  •••»  x  )  with  7:  In(u)  -^In(Rn)  denote  a  real  vector 
function  and  an  interval  vector  function  respectively.  The 
corresponding  components  are  denoted  by  f ^ (x)  =  f^(x^,  ...»  xn) 
and  7^  (x)  =  F^fx^,  . ..,  xn)  respectively  with  f^:  u  -*■  F  and 
T  ;  En(u)  -*■  E(R).  The  function  f  (x)  is  called  an  interval 
extension  of  f (x)  on  u,  if 

f(x)  :=  f(x)  for  all  x  «  u; 


and  if 


f  (x)  C  F(y )  for  x  Cy  and  for  all  x,  y  t  u  .  (6) 


If  f (x)  is  a  rational  function  then  an  interval  extension  f(x) 
of  f (x)  (the  "natural”  extension)  is  obtained  by  replacing  the 
real  variable  x  by  an  interval  variable  x  and  the  real 
arithmetic  operators  in  f (x)  by  the  corresponding  interval 
arithmetic  operators. 


Let  f:  u  ■+  F n  and  let  f  c  (u) .  Let  the  (real)  n  x  n 


matrix  <J>(x,  y)  have  the  components  <p 


yv 


%v(xl'  x2 '  xv'  yv'  yv+l'  **•'  yn)  for  M'  v  =  1(1)n'  which 

are  defined  as  divided  differences  by  (f  :=  3f  /3x  ) 

•*  yv  y  v 

r 

fy (Xj, . . . >xv_1,yv>yv+1> . . . >yn)-fu (xlf . . . >xv_1,xv,yv+1. 


•'V 


yv  -  xv 


.. 


yv 


f y v  (xl'*-*'xv-l,xv'yv+l,***,yn) 


for  xv  ?  yv 


for  xv  =  yv 
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The  matrix  <Hx,  y)  is  continuous  for  x,  y  *  u.  Define  the  class 
of  functions 

3  :=  {f(x)  t  c\u)  j  1  (x,  y)  exists  with 

y)  «  B(x,  x)  C  B* 

for  all  x,  y  t  x  C  u,  where  B(x,  x) 
has  the  property  (6)  with  respect  to  the 
variable  x  and 

B:  u  x  ln(u)  •>  In(Rn)  x  In(]Rn)  and 

B*  t  IIn(Rn)  x  En(  Rn)  }. 

3.  Results. 

It  is  in  general  difficult  to  decide  if  <j>  ^(x,  y)  exists 
for  a  given  function  f (x)  and  if  there  are  bounds  B(x,  x)  and  B* 
such  that  (7)  holds-  The  following  criterion  provides  a 
sufficient  condition  which  can  easily  be  checked  on  a  computer: 

Criterion:  Let  f  for  u,  v  =  l{l)n  be  interval  extensions  of 

f  :=  8f  /3x  on  u.  Let  the  interval  matrix 
_UXv  H-  v 

F ' (x ,  x)  have  the  components 

fyu  (xr  . xv-1,  xv,...,xn).  Let  F'  (u,  u)  be 
nonsingular.  Then  f  t  3  with  B(x,  x)  := 

F,_1(x,  x)  and  B*  :  =  F'"1 (u,  u) . 

Proof:  By  the  (one  dimensional)  mean  value  theorem  there  exist 

real  numbers  t  [min(xv,  yy) ,  max(xv,  yv)]  C  u  such  that  for 

Vt  v  =  1  (l)n 


(7) 
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’’yv  ”  ’’yv 

t  • 

•  * •  xv'  yv'  • • • '  yn} 

(x1# 

**•'  xv-l'  zy v ’  yv+l' 

yn> 

u 

t 

•••»  xv_i»  zyv*  yv+i'  •••' 

yn> 

« 

(X^ 

xv_i' vxw+1/*--.xn) 

C  f  (x)Cf  (u  )  . 

|JLV'  HI/ 

Hence  <J>(x,  y)  <  F1  (x,  x)  C  F'  (x,  x)  C  F*  (u,  u)  .  Since  F'  *  (u,  u) 
exists,  there  exist  also  the  other  inverses  with 

♦-1(xf  y)  <  B(x,  x)  :  =  F'-1  (x,  x)  CF,-1(x,  x)  C  B* 

:  =  F '  (u,  u)  . 

The  basic  idea  of  this  criterion  comes  from  a  paper  from  E.  R. 

Hansen  [2] . 

Lemma:  Let  f  «  3  and  let  x,  y  «  u.  Then  the  following  identity 

holds : 

f(y)  =  f(x)  +  <Mx,  y)  (y  -  x)  .  (8) 

The  proof  is  very  simple:  By  using  the  definition  of  the 
one  gets  for  y,  v  =  1(1) n: 

fy(y)  =  fy(yl,y2,y3,,*’,yn> 

"  fy(xl'y2'y3 . yn)  +  (yl~xl)<pyl(xl'yl'y2 . yn> 

"  f  y  (xl'x2  ,y3 '  *  *  *  ,yn)+  (yl"Xl)  *yl  (xl'yl'y2  '  *  *  *  ,yn} 

+(y2-x2)Vy2(x1,x2,y2,...,yn) 
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n 

f'V  ^  +v^  ^yv_xv^  ’’uv  *xl' —  'xv'yv'-*-'yn)  • 

These  are  the  components  of  (8). 

Definition  (Interval  Newton  Operator) : 

Let  f«  3  and  x  c  x  C  u.  Let  the  interval  Newton  operator 
N(x,  x)  with  N:  u  x  En(Rn)  +  In  (]R  n)  be  defined  by 

N(x,  x)  :=  x  -  B(x,  x)  f(x)  .  (9) 

Theorem;  Let  x  £  x  C  u,  let  f  £  g  ,  let  N(x,  x)  be  defined  by  (9) 
and  let  x  «  u  with  f (x)  =  0.  Then  the  following  5  statements 
hold: 


a)  Any  solution  x  of  (1)  in  u  is  unique  . 

b)  If  x  t  x  then  x  £  N(x,  x)  . 

c)  If  x  D  N  (x ,  x)  =  <t>,  then  x  {  x. 

d)  If  N (x ,  x)  c  x  then  x  £  x . 

e)  Let  the  real  vector  sequence  {x^}  with  x^  £  u  for 

v  =  0 ,  1 ,  ...  have  the  property 

£im  x  ^  =  x  . 
v><» 

Then  the  interval  vector  sequence  {x^},  defined  for 
v  =  0 ,  1 ,  ...  by 
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has  the  two  properties 

x  <  x(v)  for  v  =  0,  1,  ...»  (11) 

and 

£im  x(v)  =  x  .  (12) 

v-*°° 

Proof :  Ad  a)  Let  x,  y  e  u  be  two  zeros  of  f (x) .  Then  from 
(8)  it  follows  that  4>(x,  y)  (y  -  x)  =  0,  but  since  $  *  exists  for 
all  x,  y  (  u  this  implies  x  =  y . 

Ad  b)  Put  y  :=  x  in  (8).  This  gives 
0  =  f  (x)  +  <j>(x,  x)  (x  -  x)  .  Since  $-1(x,  x)  exists  and  with  (7) 
one  gets  immediately  x  =  x  -  <t>  ^  ( x ,  x)  f  (x)  t  x  -  B(xf  x)  f(x)  = 

N(x,  x)  . 

Ad  c)  Suppose  x  t  x.  Then  from  b)  it  follows  that  also 
x  «  N(x,  x)  which  contradicts  the  assumption  x  H  N(x,  x)  =  <{>. 

Hence  x  ^  x. 

Ad  d)  Define  for  fixed  x  e  x  the  real  operator 
Ty  :=  x  -  ^(x,  y)  f(x)=y  -<})"\x,y).  Clearly  T  is  continuous  for  all 

y €  x.  Furthermore  because  of  Ty  t  x  -  B (x,  x)  f (x) 

=  N(x,  x)  C  x,  the  jperator  T  maps  x  into  x.  Since  x  is  a  convex 
set,  the  Schauder  fixpoint  theorem  can  be  applied  with  the 
result:  there  exists  a  real  vector  x  c  x  such  that  x  =  Tx. 

This  implies  4>-^(x,  x)  f(x)  =  0  and  therefore  f  (x)  =  0,  i.e.  x 
is  also  a  zero  of  f (x) . 

Ad  e)  The  property  (11)  follows  immediately  from  b) . 
Furthermore  since  x^*  f  y*v^  C  u  and  B(x,  x)  has  the  property  ( 6) 
in  x  it  follows  that 
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u)  f(xlv)) 


N(xtv) ,  y tv)  >  t  N(x(v) ,  u) 

-  x<*>  -  B(x1'", 
Cx1'”  -  B*  £(xlv>)  , 


-  (y]  _  (\j)  ^ 

hence  Jiim  N(x  ,  y  )  =  x.  Therefore  with 


\)->CO 


x(v-l)  n  N(X(V)  ,  y^)  Cn(x^,  y^)  it  follows  also  that 
dim  x^v  ^  H  N(x^,  y^)  =  x,  which  completes  the  proof  of  the 

v-+-°° 

theorem. 


Remarks :  Ad  a)  The  class  3  is  so  small  that  it  contains  only 

functions  f (x)  with  at  most  one  zero  x  in  u. 

Ad  b)  Any  solution  x  of  (1)  contained  in  the  interval 
vector  x  is  also  contained  in  the  image  of  the  Newton  trans¬ 
formation:  i.e.,  the  Newton  transformation  "does  not  lose" 

any  zero  of  f (x) . 

Ad  c)  This  is  a  criterion  for  the  non-existence  of  a 
zero  x  of  f (x)  in  the  interval  vector  x.  This  criterion  can  be 
programmed  to  be  checked  by  a  computer! 

Ad  d)  This  is  a  criterion  for  the  existence  of  a 
zero  x  of  f (x)  in  the  interval  vector  x.  This  criterion  can  also 
be  checked  for  by  a  computer . 

Ad  e)  If  somehow  a  convergent  real  vector  sequence 
{x^v^}  with  x ^  x  for  v  -*■  00  has  been  determined,  the 

application  of  the  Newton  operator  by  (10)  produces  immediately  a 
sequence  (x^  }  of  bounds  to  {x^}  and  x  (a  posteriori  error 
bounds),  such  that  x^  {  x^  for  v  =  0,  1,  ...  .  This 
sequence  converges  too  toward  x.  The  speed  of  convergence 


-12- 


#1136 


(linear,  superlinear,  quadratic,  ...)  can  easily  be  related  to 
the  speed  of  convergence  of  {x  ^  . 

If  x  t  x  then  the  Newton  Operator  N  can  be  applied  to  the 
pair  x,  x  and  the  result  b)  of  the  theorem  holds.  Assume  that 
the  interval  vector  x^v~^  has  already  been  computed.  Since  the 
sequence  {x^v*}  is  defined  independently  of  the  sequence  {x^}, 
no  guarantee  is  given  that  the  next  term  x^  t  x^v  ^  and  that 
therefore  xc  N(x^,  x^v-1^).  This  is  the  reason  why  the 
intermediate  interval  vector  y^  is  introduced  in  (10). 

4 .  Example .  Convex  programming . 

Let  ip  =  \p  (x)  with  <J;:  R  n  -*■  Rand  i|>  t  C2<Rn)  and  let  ip  be 
bounded  from  below  on  Rn.  The  unconstrained  minimization 
problem 


ip(x)  -  minimum  (13) 

is  considered.  Define  the  function 
f  (x)  :=  d\p/dx 

with  f:  R  n -*•  Rn  and  f  *  C^OR11)  and  define  the  set 
S  :=  (x  t  Rn|\Mx)  <  ^(x^0^)} 

/  A  \  —  /\ 

for  a  given  real  vector  x '  «  R  .  Let  S  be  a  convex  set 

containing  S.  Let  the  Hessian  matrix 

H(x)  :«  32<J>/3x2  -  3f/3x 

satisfy  the  condition 
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W(x,  x)  <  (x,  H(y)x)  <  X (x,  x) 

with  p  >  0  (where  is  the  inner  product  in  Rn)for  any 

n  A 

vector  x  «  R  amd  y  t  S.  Then  S  is  bounded  and  the  minimum 
problem  (13)  is  in  S  a  convex  problem  and  has  exactly  one 
solution  which  is  the  only  zero  x  of  f (x)  in  S.  Furthermore 
the  real  Newton  method  (2),  (3)  applied  to  the  function  f(x) 
and  starting  with  x^  produces  a  sequence  {x^}  with  *  S 

which  always  converges  to  x  (where  perhaps  a  finite  number  of 
steps  have  to  be  changed  in  order  to  get  always  t|>(x^v+^)  < 
<Mx^)),  see  Goldstein  [1]. 

Hence  the  above  suggestion  cam  be  used,  provided  f  «  3  for 
a  certain  interval  vector  u  containing  S.  This  being  true  the 
interval  vector  sequence  {x^}  defined  by  (10)  has  the 
properties  (11)  and  (12),  i.e.,  provides  error  bounds  and 
converges  toward  x.  To  the  author*  s  knowledge  to  the  problem 
(13)  no  other  procedure  is  yet  known  which  gives  an  interval 
vector  sequence  with  the  properties  (11)  and  (12). 


Numerical  example. 

The  following  example  is  chosen  to  be  as  simple  as  possible 

in  order  to  show  all  the  steps  of  the  above  method.  Let  n  =  2, 

2  2  2  2 

t  >  0  and  2^(x^,  x2)  :  =  xi  +  x2  +  e  x^  x2.  Then 

I  xx(l  +  e  x2)' 

f(x1#  x2)  ;=  9^/ax  := 

'  x2(l  +  e  xj)/ 

and 
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3f/3x 


32<Ji/3x2 


1  +  e  xi 


2  e  xi  x2 


\ 


2  e  Xi  x2 


1  +  e  x 


1  / 


It  is  trivial  that  x^  *  x2  =  0  is  the  unique  solution  of  (13)  in 

2  _  _ 

R  .  Let  y  :=  [-1,  1]  and  let  u  have  the  components 


:=  Y- 


Define  S  :  =  {x 


ll>' 


I  x2  I  -  } .  Then 


,<°)  .=  x(0>  .= 

“1  *  x2 


u  C  s  :=  {x  f  R  2  |iji(x)  £  <|>(x*°* )  }  C  S.  One  sees  easily  that 


:=  1  >  0  and  X  :=  l+6e+3e  can  be  chosen. 


/  1  + 


( 3f/3x) 


'1  .= 


£  X. 


-2  e  Xj^  x2 


Furthermore 

\ 


y  -2  e  x1  x2  1  +  e  x2 


/det  , 


?  2  2  2  2 

where  det  :=  1  +  e (x£  +  x2)  -  3  e  x^  x2  >  0  for 

2  2  — 

max(x^,  x2)  <  1/e,  thus  e  <  1  in  u.  From  this  result  it  is 
simply  shown  that  the  real  Newton  method  (2) ,  (3) ,  starting  with 
x ^  gives  a  convergent  sequence  {x^}  with  t  u  and 


(v) 


-*•  x  for  v  -*■  <*>,  where  x 


(v+!)  :=  x(v+l)  J=  2  £  x(v)3/ 


(1  +  3  e  x|v)2)  for  v 


0/  If  •••  • 

To  apply  the  interval  Newton  method  (4),  (5),  compute 

e  x1(x2  +  y 2 ) \ 


/ 1  +  e  y2 


4>(x,  y)  :■ 


e  y2<xi  +  yj.) 


1  +  e  x 


/ 


and 
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-e  x1(x2  +  y2) 


.-1 


(x,  y)  := 


1  +  e  xi 


y2<x!  +  1  + 


e  y. 


/ 


det  <}>  (x,y ) , 


with  det  4>(x,  y)  :=  1  +  e  (x^  +  y2)  -  £  x1  y2  (Xj.  x2  +  y2  x2  +  y2* 


-e  x1(x2  +  x2)\ 


/ 


det  ij>(x,x)  , 


—  o  o  2  —  _  _  — . 

with  det  5(x,  x)  :=  1  +  e (x^  +  JJ)  ■  E  Xj  x2 ^xi  x2  +  X1  x2  +  X1  x2^' 


Finally  B(xf  x)  e  B*  :=  B(u,  u)  with 
/  1  +  e  72  2  e  y2  \ 


B* 


2  e  y2 


1  +  e  y2 


/det  B*  ,  and 


det  B*  :=  1  +  2  e  y2  +  3  e2  y4  »  1  +  e(2  +  3t)  [-1,  1]  > 


0  for 


0  ^  e  <  1/3.  Kindly  note  that  the  interval  matrix  B*  is  only 
defined  for  0  ^  e  <  1/3,  while  was  defined  in  the  much  larger 
region  0  <  e  <  1.  This  loss  is  quite  typical  for  the  switch 
from  real  to  interval  arithmetic  in  the  case  of  matrices. 

To  apply  the  above  criterion  one  computes 


/  1  +  e  x\ 


F'  (x,  x)  :> 


2  e  X1  x2 


2  e  x^  x2  1  +  e 
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Hence 


F'  1(x,  x)  := 


1  +  e  x| 


-2  e  x^  x2 


with  <3et  F  (x,  x)  (1  +  e  x^)  (1  +  e  x^)  -  4  x1  x^  ^2 

This  gives 


-2  e  x1  x2 


1  +  e  x22 


cfet  F(x,x) , 


/ 


-  -2 
x. 


F ' -1  (u,  u)  :  = 


i  x  “2 

1  +  £  Y 


o  -2 

^2  e  Y 


2  e  y2  ^ 


1  +  E  Y2/ 


det  F (u,u) , 


with  3et  F(u,  u)  =  1  +  e{2  +  5e)  [-1  1]  >  0  for  0  <  e  <  (/5-l)/5, 
By  comparing  these  results  one  sees  clearly,  as  stated  in  the 
above  criterion,  that  B(x,  x)  C  F,-1(x,x)  and  B*  C  F,-1(u,  u)  . 


Kindly  note  that  B*  is  defined  for  0  <_  t  <  1/3  while 
F '  1  (u,  u)  is  defined  only  : 

0  <  E  <  (/5-l)/5  *  0.29(1)  . 


F'  1(u,  u)  is  defined  only  in  the  still  smaller  interval 


Table  1  shows  some  numerical  results  for  e  =  0.5  and 
v  =  0(1)5.  We  remark  that  the  symmetry  of  the  problem  which  is 
illustrated  in  the  real  sequence  {x^}  by  x|v^  =  x2v^  is  not 
carried  over  to  the  sequence  {x^},  due  to  the  asymmetry  of 
B (x ,  x)  and  therefore  of  N(x,  x) .  Please  note  also,  that  for  the 
higher  values  of  v  it  is  not  true,  that  x^  <  x^l  This  means: 
the  error  interval  x^  containing  the  solution  x  is  not  only  an 
error  bound  for  x  but  also  an  improvement  over  the  approximation 
x^v\  This  very  favorable  property  is  quite  typical  for  the  use 
of  the  Newton  operator,  at  least  in  the  neighbourhood  of  the  solution. 
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Table 


)  rounded  to  three  significant  digits. 
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